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Abstract. - We study the extra term of particle current in a 2D k-cubic Rashba spin-orbit cou- 
pling system and the integer quantization of Hall conductance in this system. We provide a correct 
formula of charge current in this system and the careful consideration of extra currents provides 
a stronger theoretical basis for the theory of quantum Hall effect which has not been considered 
before. The nontrivial extra contribution to particle current density and local conductivity, which 
originates from cubic dependence on the momentum operator in the Hamiltonian, will have no 
effect on the integer quantization of Hall conductance. The extension of Noether's theorem for 
the 2D k-cubic Rashba system is also addressed. Two methods reach to exactly the same results. 
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Introduction. — Recent experimental demonstra- 
tions of spin Hall effect in some semiconductors [1-5] 
may create a way to manipulate the spin of carriers in 
terms of electric field that presents potential in future 
applications. It has stimulated many scientists' interest. 
The experiments clearly show that the spin-orbit cou- 
pling (SOC) of carriers in some semiconductors plays a 
key role in disclosing the spin Hall effect. Several theo- 
retical models of SOC have been suggested to study the 
charge and spin transport for different kinds of semicon- 
ductor systems, such as 2 dimensional (2D) linear fc depen- 
dent Rashba [6] and Dresselhaus [7, 8] models, quadratic 
fc dependent Luttinger model [9], 3D fc-cubic Dresselhaus 
model [10] and 2D fc-cubic Rashba SOC model which was 
found in a GaAs-Al^Gai-^As interface of a typical semi- 
conductor heterojunction where the fc-cubic Rashba SOC 
effect for heavy holes can not be neglected in a high- 
density regime [11]. For some Hamiltonians including 
terms with high order power (> 2) of momentum opera- 
tors (MO) p, like 3D fc-cubic Dresselhaus model, we have 
proved the conventional expression of particle current den- 
sity (CD) j corMJ (r,^ = Re{^(r, t)l/(ih)[r, H]^{v, t)} = 
(1/e) Re{^{r,t)(dH/dA)i)(v,t)} is no longer valid [12]. 
In that system, for the sake of current conserva- 
tion, a nontrivial extra term of CD j e xtra(j, t) (V- 
jextra(r, t) ^ 0) should be added to the conventional one. 
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Then the continuity equation of conserved particle CD 

ex I ra 

(r,i)) can be satisfied. Thus, 
the extra term j e xtra{^, t) is a physical quantity and has 
effect on the conductivity of the system. It is natural to 
extend it to fc-cubic Rashba system where the extra terms 
of CD may also appear due to cubic fc, the term with high 
order power of MO p in its Hamiltonian. However, 2D 
fc-cubic Rashba is a real system that can demonstrate the 
integer quantum Hall conductance. The high precision of 
integer quantum Hall effect [13] was explained in some fa- 
mous papers [14-16] where the expression of charge CD 
is implicitly based upon the conventional form. One nat- 
urally questions whether some correction to the quantum 
Hall conductance could come from the additional term of 
current j e xtra,(r, t) in 2D fc-cubic Rashba SOC semicon- 
ductors. In this paper, we would rigorously deduce the 
exact expression of charge CD that shows the existence of 
nontrivial extra term j extra (r, t) (V ■ j e xtra(r,t) ^ 0). It 
is not a local circular current and does have the contribu- 
tion to electric conductivity. Further, we prove that it has 
no contribution to the quantum Hall conductance. So the 
explanation of integer quantum Hall conductance is ex- 
tended to a more general case that includes MO of triple 
power in the Hamiltonian, though whose formula of charge 
CD must be corrected by a nontrivial extra term due to 
the requirement of its continuity. Our paper shows a more 
clear understanding of the property of integer quantization 
of Hall conductance no matter the Hamiltonian including 
additional cubic fc dependent SOC which is a realizable 
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2D quantum Hall system. 

This paper is organized as following. Firstly, we sim- 
ply introduce the formulae of the calculation of particle 
CD in the first section. The deduction of the expres- 
sion of extra term j e:C i ra (r, t) in addition to j con v(j,t) for 
a 2D cubic Rashba Hamiltonian is presented. The sec- 
ond section gives a proof that there is no contribution to 
the integer quantum Hall conductance from extra term 
jextra{r,t). Our expression of particle CD confirmed by 
extended Noether's theorem is attached in the appendix. 

Density of particle current. — We study the 2D 
cubic Rashba system that is a promising model system for 
an ultra thin film of p-doped semiconductor [11]. In a per- 
pendicular magnetic field, the single particle Hamiltonian 



H = H N (p,r)+H R , (1) 
H N (p,r) = f/2m* + V(r)-eyE y , (2) 

(3) 



H R = i\(p_a — p + <7 



where V{r) is a local spin independent potential and could 
contain an impurity potential, A — a/2h 3 is the spin-orbit 
coupling constant, E y is the transverse Hall electric field, 
p ± — p x ± ip , a = a x ± i(7 y where o x and a y are Pauli 

matrices, and p = (p x — eA x ,p y — eA y ) . The correspond- 
ing Schrodinger (or say Pauli) equation is 



(4) 



where Hamiltonian H is a 2 x 2 matrix. The particle 
density for a pure quantum state is n(r, t) = tp^(r, t)tp(r, t) 
in which we have performed the inner product for spin 
space, but not for position. This rule of inner product 
is also used in the following deductions implicitly. Since 
the number of particles is conserved, the total number 
of particles N = J n(r, t)dr should be a constant. The 
conserved particle CD j(r, t) is defined by the following- 
continuity equation: 



<9n(r, t) 

at 



-V-j(r,t). 



(5) 



For simplifying the notations, in the paper, we will not 
discriminate the notions of particle CD and charge CD 
which only differ by a factor of charge e and can be self- 
explanatory according to the context. 

For a mixed state, the density matrix p = 
^2 n lV'n}Pn {^n\ where p n is the probability of the state 
I'tpn), Pn > 0, Yl Pn = 1-The density of particle is defined 

n 

by 

n(r,t) = (r,t\p\r,t) =^2(r,t\ip n ) p n (ip n \r,t) 

n 

= ^2pi^n(r,t)tp n (r,t), (6) 



We discuss the case of p n being time independent. Then 
based on the Schrodinger equation, the left hand side of 
eq.(5) can be expressed as 



dn(r, t) 
dt 



+{-^H R ip n {r,t))^ n (r,t) 
in 

+^(r,t)(l# R Vn(r,i))}, 
in 



(7) 



where j%(r,t) = Re |W(r, t) (l/(ifi)[r, H N ]ip n {r, t)) } 
is just the conventional expression of parti- 
cle CD for non-SOC part H N (p,r). If the 
last two terms can be changed into —V • 
Re^J(r,t)(l/(iH)[r,H R ]ij n (r,t))\, the particle CD 

would bej"(r,t) = Re{i^{r,t){l/{ih)[r 7 H}4' n (r,t))}, 
which is just the conventional formula Jconv( r ^)- How- 
ever, it is not right. There should be an extra term 
jextra{ r ,t) m addition to the term i™ onv {Y,t). And we 
will prove that V • j e xtra( r , t) ^ 0. Thus, the conventional 
formula of particle CD is not conserved in 2D cubic 
Rashba system. For simplicity, we only consider the case 
in pure state and denote ip = ip n (r,t). The expression of 
particle CD for mixed state can be easily obtained from 
the one of pure state. After some algebra, we obtain 



^[r,H R ]ip\ V + ^ (^_[r,H R ]ij; 



1 eBy 
-V-jfl(r,t). 



( x,H R ip) ip + ip^ 



x, H R 



where 



-XRe 



-2A7T 



{d x 4>y <Jyd x 1p + (dytpy <Jy (dylf)) 



-i2\heBy 



+2A {eBf [ip^y 2 a y ip] 



(8) 
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t 



-[y,H R }^\ 4, + ^ l^[y,H R ]iP 



1 



-XRe[(d y ^ [y : p 2 _a + +p 2 + <j-] i>} 



XeBy , + r ~2 + ~2 -1 / 



-2A/r 



(9 x V) t o-^^V + (5 y -0) f cr x d y i> 
—i2XheBy (d y ip)^ a y tp — ip^dydyip 
~2X(eB) 2 [^y 2 a x ^\ . 
In the above, following relations are applied: 



(9) 



~2 + 

p_a 



Hr 



~2 - 

p + a 



V- [r,H F 



eBy 1 
3 ift 



x, Hi 



3 
1 

-eBy— [a;,p 2 _ < T + +p^cr-] . 
2 ift L J 



So we have j(r, t) 
term is 



= jjv(r,i) + jfl(r, i). Then the extra 



3extra{r,t) — j (r, £J — jconu (r, t) 

= j R (v,t)-Re^^[r,H R ]iP 



where j c 



,t) = i?e{'0 + l/(i^)[r, i/]^} is the so called 



conventional current that widely appeared in literatures. 
Then the expression of extra particle CD can be finally 
simplified as 

f extra (r, t) = 2Xh 2 d x d y tyW) _ (^a y i>) 

+Xh 2 d 2 y , (10) 

f extra (r, t) = 2Xh 2 d x d y tyV^) + Xh 2 dl 

-Xh 2 d 2 . (11) 

In above equations, all (f/'Vo.V') = ( r j t), a x tp(r,t)) 
are position dependent. The divergence of 

jexi™(r, t) is generally non-zero, V ■ j extr a(*,t) = 
-XK 2 dl (^a y ^) - XH 2 dy (^a x ^)+ih 2 Xd x d 2 y (i(>* a y i/>) + 
'SXh 2 d 2 d y Ctp^a x ip) ^ 0. So, as shown in Eqs.(10) and 
(11), we derived non-trivial extra terms in the expression 
of conserved particle CD of a 2D cubic Rashba Hamil- 
tonian. The same result of extra currents can also be 
obtained by taking account of the gauge invariance based 
on Noethcr's theorem. Its detail is presented in appendix. 

When we consider a mixed state, the extra term of 
charge CD can be expressed as 



iextra( r ,t) — e ^ ] PnjeJtra ( r ; ^) • 



(12) 



In fact, our Hamiltonian H{= Hn + H R ) is time inde- 
pendent, so the particle density n and charge CD j are 
position dependent only. 



Hall conductance. — Now we study the charge CD 
along x direction which is j (x) (r) — jconv (r) + j extra ( r ) ■ 
We take the integral with respect to y for j^ x \r) to get 
the charge current 

I (x) =! j£U*)dy+ [ j£L(r)dy, (13) 



where L y is the width of the system. And denote the 
length of the system as L x . Finally L x and L y can ap- 
proach to infinity if the system becomes macroscopic. I^ x ' 
should not be position x dependent because of the particle 
conservation. Then we take the integral of x for I^ x ' : 



jO) = _ 



l ~ [ [ jtLWdr+i f f jtL(r)d: 



t{x) , t(x) 
conv * extra' 



where = L x L y , I^lv = ^ JJ jcoL{r)L y dr and tf x \ ra 



(14) 



(x) 



r(*) 



H // jextra( T )Lydr. From Eq.(10), the extra part of the 



current is 



r(*) 



eXh 2 L y J2Pn{^ f f [2d x dy (*t(r)o- x * n (r)) 
a 

-dl (*t (r)a y * n (r)) + d 2 y (*t ( r )<r s * n (r))]dr} 

-eXLy^2p n ^ JJ dY{2{p x p y ^ n {r))^ a x ^ n (r) 

n n 

+2*t (r)a x p x p y ^ n (r) - 2 (p x ij' ri (r)) t a x p y ^ n (r) 
-2(p y *„(r)) t a x p x y n (r) - (p^ n (r)) f a y ^ n {r) 
-*3»(r)o-»^*n(r) + 2 (^*„(r)) f <r v p B *„(r) 
+ (pliWj'^iW + ^ n (r)<Jyp 2 y ^ n (v) 
-2(p y ^ n (r)^a y py^ n (r)}. (15) 



The terms in right side of the above equation become the 
spatial inner product after the integration of r over the 
whole space of the system. Since the operators {p x ,p y } 
are hermitian, as an example, we have 

^JJ dr{2(p x py^ n {r)) ] <7 x * n (r)} 
n 

= i JJ dr{2(p v ^ n {r)) f a x p x ^ n (r)} 
n 

= ^ JJ dr{2 (*„(r)) f ^^^(r)}. 
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Considering above property for inner product of position 
space in equation (15), we can easily obtain 7 ( 



0*0 



0. 



No contribution to Hall conductance from extra term of 
charge CD is proved. Finally, wc have 



Ax) 

Jconv 



(r)dr 



n n 
eL y ReJ2pn (*«| ^z[x,H ] |* n ) , 

n 

p 2 /(2m* 



i\(p_a + — p + o ) 



Thus, the quantum Hall conductance is only from the con- 
ventional term j co ™(r, t). For Ho, the cubic 2D Rashba 
model without transverse electric field, its Schrodinger 
equation is 

H \^)=E^\^). 
It has been solved exactly [17], 

£4 0) = (n+l/2)hw,n<2, 

E { J°l = (n-l) + s\h 2 n(n-l)(n-2) + - fnu, 



(16) 



(n - 1) + sy 7 2 «(« - - 2) + - 
71 > 3. 



where u> = eB / (toc), s = ±1 and 7 = 4Ato* ^/2heB/c. 
The cigen energies in traditional quantum Hall effect are 
Landau levels separated by gaps. Now the "Landau levels" 
of a 2D cubic Rashba model have some modification for 
?i > 3, but they still keep the essential feature of the gap 
separation. The corresponding cigenfunctions are 



l*£ 0) > 
l^ 0) > 







, n < 2, 



1*5) 



C n sl4>n-3, 
C ns 2(f>n 



,n > 3. (17) 



where {C ns \, C nS 2] are normalized constants, 



s2 



1 



"fy/n(n — l)(n — 2) 



J + S J 7 2 n(n-l)(n-2) + || 



and is the wave function of harmonic oscillation type. 
Impurities may result in widening out the "Landau lev- 
els" . The conventional velocity operator in position space 
is v = 1/(2/1) [r, H] = l/{ih) [r, Hq]. It is easy to have 
v(k) = l//W kJ E (0) (k) in k space. Following Laughlin [14] 
or Kohmoto's [16] deduction, the integer quantization of 



quantum Hall conductance can be obtained. Here we will 
present a different approach to reach to the conclusion of 
integer quantum Hall conductance for such a specific 2D 
cubic Rashba system. 

Since the cigen energies and wave functions of 
Schrodinger equation in the second quantization represen- 
tation can be found exactly, we also calculate the Hall 
conductance in linear response approximation and it shows 
excellent consistency with integer quantization of Hall con- 
ductance. More specifically, the Hall conductance a xy of 
this system can be written as a xy = J2 n .s N n,s ( a xy) ns , 
where N„ tS is the number of particles occupying the (n, s)- 
th "Landau level" (here the "Landau level" is marked by 
two index, n indicating the energy level of the system with- 
out SOC, s indicating the energy level splitting due to 
SOC). And (&xy) n s is the one particle's contribution from 
the (71, s)-th "Landau level" , by linear response theory, 



( a xv) n ,s ~ ^ 

(n" ,s")^(n,s) 





jx 






H' 




EyL X Ly 




.,) 





h.c{\lS) 



where H' = —eyE y and the electric field is uniform. Now 
we adopt the Landau gauge p x — hk x — eBy/c,p x ~ p y , 
and introduce the operator of bosonic quasi particles a = 

(fx - iPy) , flt = VlS (p* + l Py) ' [«' a 1 = L 



V 7 



2heB 

Then we get 



H Q = flUJ 



a' a + 



17 a 



Then 



Jx 



e 
ih 



x,H 



he , +x 3ie7 heB 

(a + gJ) H -\ 

2eB K ' to V 2c 



it 2 



And 



W 



-eyEy 



-eE y hk x 



he 
2e~B 



(a + ot) (20) 



where / is a unit matrix. Then using Eqs.(17), (19) 



and (20) for matrix elements ilisl^^+i ±) and 



v(0) 



(^i il^'l^n'+i ±)i ^ ne selection rules will be found. And 
the summation over states in Hall conductance in Eq.(18) 



,(()) 



T1_/l 



Extra current and IQH conductance in the SOC system 



can be simplified as 



1 



(0) 
n,± 



(0) 

n+l,+ 



(0) 

n+l,+ 



(0) 



EL x L y 



E 



(0) 



n,± 



3x 



(0) 



(0) 

n—l,+ 



,(0) 
'n+1, 

w 



(0) 
n,± 



(*i 0) ± 




n+l, — 1 \ x n+1,— 




*&) 




jx 


P (0) p (0) 

$(°), \/* ( °\ 

n — 1 , — / \ n — 1, — 




*i 0) +) 



(0) 
n,± 



Then 



(0) 



using 
3 



eqs.(16) and 
A A and /*W 



(o) 

n-l,- 

(17), 



all 



(0) 



+ h.c] 



the elements 
can be cal- 



n±l,± / clliu \ *n±l,± 
culated without difficulty. After long but straight alge 
braic deduction, we can finally obtain 



(<7xy) 



(21) 



where <£>o = BL x L y ,$ = hc/e. By summing up all the 
contributions from different energy levels, the total Hall 
conductance will be 

e 2 1 



n n.s 

e 2 N Q 



h <&/$o V h 



h $/$ 
(22) 



Here N is the total number of carriers, and the filling 
factor v = AT /($/$ ). 

Due to the existence of impurities in practical samples, 
localized states appear in the region between "Landau lev- 
els" . It leads to the appearance of the plateaus when the 
Fermi level lies in that region. The gap between two con- 
ductance plateaus is obviously e 2 /h. It is concluded that 
the cubic SOC do induce the extra term of CCD that yields 
the contribution to electric conductivity, but no contribu- 
tion is given to the quantum Hall conductance. 

Conclusions. — We have derived an exact formula 
of particle current density for a 2D cubic Rashba model 
that appears in some p-doped semiconductors. In addition 
to the conventional current expression, there must be an 
extra term that ensures the current continuity equation. 
The extra term must have the contribution to electric con- 
ductivity, but no contribution to the charge quantum Hall 
conductance that is proved rigorously. So, it can be clearly 
shown that no effect is made on the topological property 
of integer quantization of Hall conductance due to the ex- 
istence of extra terms in the 2D cubic Rashba coupling 
system. Further experimentally detectable effects of the 
new term arc still on research. 



Appendix: Deduction of extra terms from 
Noether's theorem. — In this appendix, we point 
out that, for a 2D cubic Rashba Hamiltonian where the 
highest order of derivatives is higher than 2, it is neces- 
sary to generalize the expression of conserved current in 
Noether's theorem. Applying the generalized Noether's 
theorem [12], we can get the expressions of conserved par- 
ticle CD of a /c-cubic Rashba SOC system from U (1) gauge 
invariance. 

Noether's theorem, not only indicates the relation be- 
tween conserved currents and symmetries of Lagrangian, 
but also implies that the expression of conserved current 
depends on the form of Lagrangian from the begin- 
ning of its deduction. In usual cases, Lagrangians are 
expressed as C[<f) (x) , d^cj) ( x ) , $ ( x ) > 9^ (re)], x 11 = 
(t,r),/j, = 0,1,2,3 -such as the Lagrangian of com- 
plex scalar field- which only include fields <f> (x) , <jy (x) 
and their first order derivatives d^<p (x) , 8^(jr (x) 
as independent variables. But in our case, Hamil- 
tonian H includes higher order derivatives. So 
its Lagrangian should be written in the form 
C[(j>(x), d^(x), dftdv^x), ft (x), d fl ft(x),d t _ l 8„ft(x), ...],| 
where higher-order derivatives are also included 
as independent variables. For simplicity, we de- 
note <p(x) and ft (x) as <fr and ft. The Hamilto- 
nian of a fc-cubic Rashba system studied here is 
Hji = p 2 1 (2m) + iX (p^_cr + — p+<7 _ ). The corresponding 
Lagrangian can be 

C[ft 8^, 8^ ft dfaft ft ,8^,8^,818,, ft] 

= ft W ) + _L^t(^ + a 2)^ 

+2i\ft (<7 x d 3 y + a y dl) 4> 
-2i\ft {3a x d 2 d y + 3a y 8 x 8 2 ) ft 

According to the least action principle, one can easily ob- 
tain an Euler-Lagrange equation 



= -r- - d u 



d{8^) 



dad u 



dC 



dld v - 



dC 



8(8^) " w d(d*d v <i>y 



which yields the Schrodinger equation. Actually, the first 
two terms on the right-hand side of the above equation 
give the conventional formula of particle CD. The remain- 
ing parts lead to the extra terms. The corresponding con- 



served current is 



8C 



-6(d v ft-(8„- 



d{d M <f>) T 8{8^8 v ft K ~"~ r ' y ""8(d^ft 
8C 8L 

S (8^ ft - (fl/, / Q ao 



8 {8^ft 



-01 



8C 



9(8^ ft 



'd(d^ft 

)<50 + (0 ^ 0*) 



(A.l) 



which satisfies the continuity equation d^F^ = 0. We 
concentrate on the deduction of conserved particle cur- 
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rent corresponding to U (1) gauge symmetry. From in- 
finitesimal variation of fields 5(f) = ia(f),8(fr = —ia<p', the 
expression of conserved particle CD for a fc-cubic Rashba 
system is 

\2m J \2m J 
-2A[0V (d*4>) + (9^t) ovj, 

- (d x tf) a y {d x 4>)\ + 6A[0t (o*d x 8y<f>) 

- (9 x 0t) ( a *d y ct>) + (fig^t) ( o *0)] j 

j- - -"-* , (&) + (£*y* 

-2A[0V (5g0) + (9^ f ) o*0 

- (0^) ^ (a y 0)] + 6A[^ (^a x 9^) 

Comparing the above formulae with the conventional one 
jcont) = Re{(f)^ (i [r, Hfl] <?!>)}, we get the extra term of 

particle CD Jextra = j ~ jconv'- 

Jtxtra = "A^ (^O*0) + 6A (C^t) a* (9^) 

+6A (d x d y tf) (a 9 (ft) - 6A (flg^t) ((T y0) 

+3A0t ( CT ^) + 3 A (d 2 y ^) (a^) , (A.2) 
JLt™ = -H 2 + 6A (dytf) ay (d x cp) 

+6A (5*5^) (o*0) - 6A (a 2 ^) (cr x 0) 

+3A0 t (^9 2 0) + 3A (d 2 ^) (<rV) . (A.3) 

Further, it is not difficult to check that the extra term 
jextra deduced here by extended Noether's theorem and 
jextra in the second section do satisfy the equation V • 
(jextra ~ Jextra) = 0. Thus we conclude that our result of 
extra term is rigorous. 

* * * 
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